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Abstract 

In this paper we introduce the generahzed osciUator model (GOM) 
as a family of exactly solvable models useful to investigate theoret- 
ical aspects related to the statistical description of the aging state. 
GOMs are defined by a potential function V{x) and characterized by a 
zero-temperature relaxation determined by entropy barriers and par- 
tial equilibration. Analytic expressions for the effective temperature 
can be derived using a fluctuation theorem valid in the aging regime 
without the need to solve the dynamical equations for correlations and 
responses. Two classes of models are investigated in detail: the ho- 
mogeneous potential model with V{x) = {k/2p)x'^^ {p being a positive 
integer) and the wedge potential model {V{x) — k\x\) where V{x) has 
a singularity at the ground state coordinate x = 0. For the latter, 
we present some numerical simulations that reinforce the validity of 
the main analytical results. GOMs offer a conceptual framework to 
develop a statistical description of the spontaneous relaxation process 
that has been recently proposed H to be at the root of the intermit- 
tency phenomenon observed in glasses and colloids. 

1 Introduction 

Non-equilibrium phenomena is a field of research of much current inter- 
est. From turbulence in liquids to heat convection inside stars a plethora 
of systems show a very rich and complex behavior, rarely describable in 
terms of few variables The opposite situation is encountered in equilib- 
rium systems. There few parameters are needed to characterize the equilib- 
rium state and its fluctuations. Entropy, a key concept in thermodynamics, 
admits a statistical interpretation in terms of the microscopic motion of 
molecules. Boltzmann established the bridge linking the microscopic and 
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the macroscopic worlds, the central result in his theory being the relation 
S = kB^ogiW) where W is the number of configurations available to the 
system. The extension of this approach to non-equilibrium systems and the 
characterization of their behavior in terms of a few number of parameters 
still represents a major theoretical challenge. 

Two categories of non-equilibrium systems have received considerable 
attention in past years: systems in steady states and glassy systems. The 
first category encompasses all those systems driven out of equilibrium to a 
stationary state by the action of an external perturbation. The most com- 
mon example is a wire of metal with extremes in contact with two thermal 
sources at different temperatures. In this case, and if the temperature differ- 
ence is not too large, the flow of heat from the hotter to the colder source is 
described by the Fourier law. The second category encompasses all systems 
that are not in a stationary state but which properties change very slowly 
with time. Structural glasses (such as ordinary window glass) are prototype 
examples. 

The glass state is characterized by a very slow relaxation towards equilib- 
rium and by a exceedingly low rate of the energy released from the system to 
the bath during the relaxation. A useful parameter to characterize the glass 
state is the age of the glass (also called waiting time in several experimental 
protocols) which is the time elapsed since the system was prepared in the 
non-equilibrium state. In the aging state correlation functions tend to decay 
in a timescale that is of the order of the age of the system. A statistical 
interpretation of the aging state has not yet been accomplished, however 
our present understanding might not be far from resolving several of the 
most important clues. Recent experiments have observed the existence of 
intermittent fluctuations [21 |3j that could be related to dissipative processes 
characteristic of the glass state |1J El . A thermodynamic description of such 
processes could provide an important step in that direction. 

Statistical models have been always an important source of inspira- 
tion and ideas to understand the glass transition. Many types of models 
have been considered in the past, from phenomenological two-state systems 
to spin glasses, passing through a wide range of family systems such as 
hard-spheres liquids, Lennard-Jones systems, lattice models, kinetically con- 
strained models, among others. Most of these models have been investigated 
using approximate methods or numerical simulations. 

To this list we should add exactly solvable models. From a historic 
point of view, these have played an important role in the early days of 
statistical mechanics. From the urn models introduced by the Ehrenfests 
aiming to understand concepts such as entropy and thermal equilibrium, 
to the Ising model that describes phase transitions and critical phenomena, 
solvable models offer conceptual frameworks to contrast ideas and check their 
consistency by evaluating specific predictions. In this way, exactly solvable 
models have also contributed to our current understanding of glassy systems. 
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The goal of this paper is to introduce a general family of exactly solv- 
able models that might help to better understand the mechanisms behind 
the slow relaxation observed in glassy systems. We introduce the general- 
ized oscillator model (GOM) as a generalized version of a previous model 
introduced by the author These share some properties with kinetically 
constrained models as statics is trivial but dynamics is not. Therefore they 
belong to a large category of models whose dynamics can be very rich de- 
spite of the fact that the energy landscape has trivial properties. The main 
objective all along this paper is to discuss a statistical approach to the 
aging state by emphasizing the relevance of the concept of effective temper- 
ature as a useful way to quantify violations of the fluctuation-dissipation 
theorem [3 IE] • In the framework of the GOM we stress the relevance of the 
effective temperature to quantify the spectrum of intermittent energy fluc- 
tuations in the aging regime that have been experimentally observed [21 El- 
The link between intermittent effects in non-equilibrium systems and sta- 
tistical effective temperatures has been recently proposed in the framework 
of simple models for the glass transition 0J|5]. The present paper extends 
these considerations to the GOM. 



2 The generalized oscillator model (GOM) 

Generalized oscillator models (GOMs) consist of a set of non-interacting os- 
cillators each described by a continuous variable Xi and the energy function, 

N 

E = ^Vix,) , (1) 

i=l 

where V{x) is a real valued potential energy function that diverges to -|-oo 
in the limit |x| oo. For instance, the potential can be of the type V{x) = 
^x^P with p an integer value. This is called the homogeneous potential 
model, the case p = 1 corresponding to the harmonic case introduced in [^. 
Here we will only deal with potential energy functions such that the partition 
of an individual oscillator, Zi, remains finite at finite temperatures, 

Zi = r exp(-/3y(x)) , (2) 

with /3 = l/ksT, ks being the Boltzmann constant (that we will set equal 
to one) and T is the temperature of the bath with which the system is put 
in contact. 

We consider a dynamics where all oscillators are updated in parallel 
according to the rule 

Xi > Xi H ■j= , (3) 



3 



and the rj are uncorrelated Gaussian variables with = and variance 
rjXj = A^(^jj. The updating of all oscillators is carried out in parallel in a 
single move. The move is accepted according to the Metropolis rule. We 
will focus our analysis on the zero-temperature dynamics as this is the case 
where relaxation is fully determined by entropic effects. Indeed, at T = 
activated jumps over energy barriers are suppressed and relaxation proceeds 
only through search of favorable directions in phase space where the energy 
decreases. As time goes on, dynamics slows down as most of the proposed 
moves tend to push the system away from the ground state by increasing 
its energy. Only those moves that decrease the energy are accepted. Glassy 
dynamics in the GOM is then consequence of the quick decrease with time 
of the acceptance rate. 

The GOM shares with kinetically constrained models the property 
that, while statics is trivial, dynamics is complicated. Oscillators in are 
non-interacting and therefore the model has trivial static properties, the 
A^-oscillator partition function being given by Zj^ = (Zi)^ . The thermody- 
namic properties are then derived by computing the integral (HJ. 

To solve the dynamics of the GOM we have to compute the probability 
distribution of energy changes P(AE). This is defined as the probability 
density that in a given move the energy changes by an amount AE. In 
general, this quantity shows a complicated dependence on the A^-oscillator 
probability density Vt{{xi}) that the system occupies the configuration {xi} 
at time t. However, the GOM has the good property that the P{AE) de- 
pends on a finite number of observables rather than on the full configu- 
rational probability distribution. This property is characteristic of mean- 
field systems, the GOM being just another example. This makes the model 
amenable of analytical computations and a good laboratory to test many 
results regarding the glassy regime. To compute the probability P{AE) we 
consider the change of energy in an elementary move, 

N oo . N 

^^ = E[^(-^ + ;^)-M=E^i]^E^^'H-0^' • (4) 

As dynamics is stochastic the quantity AE' is a random variable whose 
distribution can be reconstructed from the moments (AE)^. An explicit 
calculation of such moments shows that only the first two moments give a 
finite contribution in the large N limit. Therefore, in the thermodynamic 
limit, P{AE) is a Gaussian distribution 

P{AE) = {2'K(t\e)~^ exp 

with mean Mab and variance (t\^ given by, 

Mae = AE = ^V"{x) , (6) 



{AE - M^Ef 



(5) 
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alj, = iAE)^-{AEf = A'iV\xW , (7) 

where /(x) = (1/A^) J2i^i fi^i)- At T = the equations for the acceptance 
a (i.e. the fraction of accepted moves) and the energy per oscillator e = 
E/N = V{x) can be written as, 



de 
dt 



P{y)dy , (8) 

oo 

yP{y)dy . (9) 



Inserting ((5)) in H8I9|1 we obtain, 

1 . )\/r2 



at- UvrJ 2alJ^ 2 ^^^^l (2^|^)i J ' ^^"^ 



^Erfc[ 



Mae 



{2a: ' 



2 , 2 



(11) 

with Erfc{x) = {2/^/^^) ex.p{—u'^)du the complementary error function. 
These equations are not generally solvable as they are not closed and the 
time evolution for the mean Mae and the variance cr"^^ is unknown. Only 
for some types of potentials V{x) a closed solution can be found. Oscillator 
models have the interesting property that relaxation at T = gets steadily 
slower as the ground state is approached. In fact, the ground state configu- 
ration {xp^} is characterized by the fact that the energy ^ is an absolute 
minimum, therefore xp^ = x*^^, with V'{x^^) = 0. According to the 
value of o"^^ steadily decreases as the ground state is approached. From 
(|10I1H) this leads to a quick decrease of the decay rate of the energy and the 
acceptance rate. In such conditions the scenario of partial equilibration (as 
described in [S] and below in Sec. El) holds. A salient feature of (|lfl|llj) is 
the dependence of the dynamical equations upon the following parameter, 

which has the dimensions of an energy (or temperature). Indeed, at T = 
in the large-time limit, the quantity A vanishes asymptotically. If we define 
X = \/Ma_e/(4A) we can then expand the complementary error function in 
(|1()|11|) around x = oo. 

Substituting this expansion in we get for the time evolution of the 
energy, 

^ /(Tij;\| exp(-x^) 
dt V Svr / x2 ■ ^ ' 
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This equation can be asymptotically solved for a quite broad family of mod- 
els. In general, the asymptotic decay of the energy can be expressed in 
terms of the parameter A (|12|) by knowing the analytic behavior of V{x) 
in the vicinity of x = 0. We show below how the parameter A in ()12[) 
plays the role of an effective temperature that quantifies violations of the 
fluctuation-dissipation theorem [Hj. 

3 Partial equilibration 

To better understand what partial equilibration means we will consider the 
case of an harmonic well V{x) = {l/2)kx'^ where k is the stiffness constant of 
the well. This case corresponds to the linear harmonic oscillator introduced 
in |S] and studied in detail in other works |1U1 lll| . For the harmonic case 
the energy is quadratic in the variables Xi, 



Eqs. (jSEI give Mae = kA'^/2,aj^j^ = k'^A^x'^ = 2kA^E/N = 2kA'^e where 
e = E/N is the energy per oscillator. We now follow the discussion pre- 
sented in Ref. The constant energy surface can be represented by an 
hypersurface centered around the origin Xi = x*^^ = (depicted as O) of 
radius R = y/2E/k. In Figure ^ we depict a schematic representation of 
the motion of a representative configuration {x^} (depicted as P) of energy 
E in phase space. The smaller dashed circle represents the region of points 
accessible from {x-*} according to the dynamics ©. All accessible points 
{xi} satisfy J2i{^i ~ ^i*)^ = i-^- lie at a distance A from {x-*} which is 
the radius of the smaller dashed circle. The accessible configurations in a 
single move lie in a spherical hypersurface of dimension N — 2 corresponding 
to the intersection of the hypersurface of energy E' and the smaller spherical 
hypersurface of radius A. We call this region the intersecting region / as 
shown in Figure ^ The final configurations contained in / lie at a distance 
R' = yj2E' jK to the origin O. The change in energy associated to this 
transition is AE = E' — E. The probability of this jump is therefore pro- 
portional to the surface of the intersecting region, P{AE) oc C^~'^ where 
C is the radius of the intersecting region. The computation of C is quite 
straightforward as can be deduced from the triangle including the points 
P, O as vertexes and whose three sides are R,R',A. In terms of R, R' and 
A, the distance C is given by the relation: R = \/ R'"^ — + V A^ — C^. 
Expressed in terms of E, AE, A we have, 



The surface il.{E,AE) corresponding to the region / of radius C, relative 




(15) 




(16) 
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Figure 1: Geometrical construction to compute P{AE). The thick lines de- 
note the departing and final energy hypersurfaces centered at O. The dashed 
circle indicates the hypersurface accessible from point P. The intersecting 
region between the accessible hypersphere centered at P and the final hy- 
persurface of energy E' defines a hypersurface I of radius C (the radius is 
represented by a thick line). See the text for a more detailed explanation. 



to the energy E of the reference configuration {x-*}, is 
n{E,AE) oc 



8E \ k 



(17) 



Using the fact that E is extensive with N this expression can be rewritten 
as, 

(^AE-^f 



n{E, AE) (X exp 



4{E/N)KA^ 



(18) 



which is then proportional to the probability distribution ((SJ- 

This construction then provides a geometric way to determine the Gaus- 
sian distribution From (jlOj) we see that C (and therefore also ^}{E, AE) 
or P{AE)) has a maximum for AE = kA'^/2. The Gaussian distribution is 
depicted in Fig. |2I Consider now a T = dynamics where only moves with 
AE < are accepted. In this case, as the relative number of configurations 
that are accessible from P goes like C^~^, the largest number of accessible 
configurations lie in the vicinity of P. Because the radius of the small dashed 
hypersphere in Fig. ^ is equal to A but the radius R of the hypersurface of 
energy E is proportional to A^^/^, dynamics constraints the system to move 
along the constant energy hypersurface in the thermodynamic limit . The 
system has then time to diffuse throughout a given energy shell of finite 
width before leaving that shell towards lower energy surfaces. This scenario 
was called partial equilibration in 
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AE 

Figure 2: Probability distribution of energy changes for the GOM as derived 
from © and H18() . The major part of energy changes are with AE > 0. At 
T = the variance of distribution o"^^ decreases relative to the mean Mae 
(i.e. the parameter A in ((T^ decreases as the ground state is approached). 
The fraction of accepted moves (dashed area with AE < 0) steadily de- 
creases with time. 

This derivation can be generalized to the GOM where the constant 
energy hypersurface is not necessarily a sphere. The crucial point in the 
argument is then the fact that A is a finite quantity while the typical spatial 
dimensions of the constant energy hypersurface are of the order \/iV- In the 
thermodynamic limit the hypersurface is locally a sphere and the length C 
(|16l) can be mathematically expressed in terms of the mean energy curvature 
V"{x) and the modulus of the energy gradient vector given by V'{x)^. The 
final result is again the most general expressions (|5p(ip7jl . 

4 The fluctuation theorem (FT) and effective tem- 
peratures 

In the GOM a scenario of partial equilibration takes place in the constant 
energy hypersurface. This is exemplified in Fig. ^ for the harmonic case 
where it is shown how the system is constrained to dwell in the constant 
energy hypersurface before moving to lower energy configurations. In this 
case, because dynamics is microscopically reversible and ergodic along the 
constant energy hypersurface, a quasi-stationary dynamics emerges with a 
probabilistic description that can be done in terms of a microcanonical mea- 
sure. However, there are important differences between this dynamical mea- 
sure and the usual microcanonical measure for equilibrium systems. While 
in the latter the energy is kept strictly constant and there is no time depen- 
dence of the microstate probability distribution, in the former the measure 
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is dynamical and time dependent. 

A key concept in such description is the notion of configurational en- 
tropy (sometimes also called complexity) S{E) which counts the number of 
configurations with a given energy E. This quantity is precisely given by 
P{AE) except for the fact that P{AE) counts the number of configurations 
of energy E' which are accessible from a reference configuration of energy 
E. In a micro canonical description of the aging state the probability to 
visit configurations with energy E' is then proportional to the number of 
configurations with that energy Q{E'), 

P{AE) oc !^(^') = exp(s{E')) , (19) 



with AE = E' — E, E' being the final energy and E the departure energy. 
It is important to emphasize that P{AE) in ()19|) (as well as in ©) is a 
probability rather than a rate. Transition rates are transition probabilities 
per unit of time, therefore have natural dimensions of frequency. We will 
denote them by W{AE). In general we can write, 

W{AE) = , (20) 

where t{E) denotes the average time to escape from a configuration of energy 
E. As these are microscopic transitions, the elementary escape process can 
be assumed to be a Poisson process described by its characteristic time t{E). 
We can now imagine a situation where, in the large N limit, the average 
escape time is identical at both energies E, E' whenever the energy difference 
AE ~ 0(1). In this case, the timescale drops from the ratio between the 
forward {E E') and reverse {E' E) transitions, 

W{AE) _ P{AE) _ ^21) 



W{-AE) P{-AE) 

Therefore, the rates W satisfy the same micorcanonical relation between 
the forward and the reverse paths as do the transition probabilities P{AE). 
Using H19|) and expanding the term in the exponential up to the first order 
term in AE we find, 



exp 



AE 



(22) 



w 

where we have defined the time-dependent effective temperature, 

T|T(t^) V dE )E=E{t^) ■ ^ 
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The effective temperature is a quantity that depends on the age of the sys- 
tem tuj through the time-dependent value of the energy E{tw)- A remark 
is now in place. The decomposition (|20)) is reminiscent of the rates used 
in trap models 1121 • However there is an important difference between the 
GOM and trap models. In the latter the timescale t{E) also depends on the 
energy of the trap but, contrarily to the present case, the rate is modified for 
energy changes AE ~ 0{1). Moreover, the partial equilibration scenario is 
difficult to visualize due to the absence of a proper configurational space. Al- 
though phenomenological trap models are very useful models to understand 
many aspects of the aging dynamics, several issues still remain controversial 
specially regarding the physical significance of FD violations jl3l 1141 [T3] . 

Relations describing ratios between transition rates for forward and re- 
verse processes in non-equilibrium systems are commonly known as fluctua- 
tion theorems (FTs). The relations (|2H22|) show a strong resemblance with 
some of these theorems, however there are important differences that we 
want to highlight. There are two general classes of fluctuation theorems. In 
the first class there are the so called entropy production FTs in stationary 
systems where the relation between forward and reverse transitions is re- 
lated to the entropy production in the asymptotic limit of large times JHl • 
In the second class there are exact non-equilibrium work relations valid at 
all times between the forward and reverse rates whenever the system is ar- 
bitrarily perturbed away from an initial equilibrium state along both the 
forward and reverse paths \17\ I18j. The most important difference between 
H21I22|1 and these theorems concerns the fact that the aging state is neither 
in an stationary state (first class) and the system does not start from an 
initial equilibrium state (second class). Moreover, the relations (|21I22() are 
not valid for energy changes AE arbitrarily separated in time as the effective 
temperature (^5]) is age dependent. 

There is another important feature of ()21|) that must be emphasized. 
Standard FTs allow for transitions to occur in both directions and certainly 
this is what the identity 1)21(1 seems to imply. Because no work is exerted 
upon the system during the relaxation, the energy change is related to the 
heat transferred between system and bath. Therefore ()2ip relates transition 
rates between identical amounts of heat that are absorbed and released 
between the system and bath. However, we face the problem that at zero 
temperature no heat can be absorbed by the system (i.e the energy can 
never increase) and therefore (|21() cannot hold. The resolution of this issue 
concerns the true meaning of the effective temperature discussed in the next 
section. 

4.1 Spontaneous relaxation and effective temperatures 

Developments during recent years in the theory of spin glasses and glasses 
have shown that aging systems show violations of the fluctuation-dissipation 
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theorem (Hj that can be quantified in terms of an effective temperature. This 
is usually defined in terms of the fluctuation-dissipation ratio (FDR) jj^l > 



where we assume that {t — tm) /t^ ^ 0{1) . In general the effective tempera- 
ture is a quantity that depends on the measured observable and the probed 
frequency u; ~ l/(t — t^) relative to the age, ujt^ = tu]/{t — t^)- The inter- 
esting meaning of H24|) is found in the low- frequency regime ojtw < < 1 where 
violations are expected to be strong. The effective temperature defined in 
this way requires the measure of the response function, i.e. the application 
of an external perturbation or field that shifts the energy levels and exerts 
mechanical work upon the system. In this case, the exerted work might ac- 
count for part of the energy transferred from the system to the bath making 
zero-field transitions with /S.E > accessible. Therefore, H21I22() have to be 
understood as the proper way of quantifying forward and reverse transitions 
in a configurational space that has been perturbed by the action of an exter- 
nal field. This establishes a way to evaluate the effective temperature from 
the rate of heat exchange between system and bath without the explicit need 
to introduce an applied external field. Note that an external field is usually 
required to evaluate the response function R{t,tw) ^■ 

The existence of the effective temperature is then related to the presence 
of a heat exchange process between system and bath that we call sponta- 
neous as it is determined by the fact that the system has been prepared in 
a non-equilibrium state. The spontanoeus relaxation is different from the 
heat exchange process (that we call stimulated) between system and bath 
typical of equilibrium systems. In particular, the stimulated process is a 
high-frequency process characterized a Gaussian distribution of exchange 
events while the spontaneous process is a low-frequency process that mani- 
fests in the form of some tails in the heat-exhanged distribution which width 
is age dependent ^. For a recent discussion of these ideas see The exis- 

tence of these two heat exchange processes is at the root of the intermittency 
phenomenon recently observed in glasses and colloids |21 El ■ 

The effective temperature can be computed in the GOM by using p2|) 
and (ESI)- Indeed, from and © we get. 



^In another context, numerical methods to compute the response function have been 
recently proposed (201 121| . 

^The words "spontaneous" and "stimulated" make explicit reference to the problem of 
light emission by atoms in a bath of photons. In that case, the spontaneous process is the 
emission of radiation by atoms in an excited state independently of the presence of the 
bath. The stimulated process, though, is the emission and absorption of energy by atoms 
induced by the bath of photons. 




t > t 



(24) 




(25) 
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and therefore, 

Te7(tJ = A(t.) = , (26) 

V (x) 

where A is an age dependent quantity that has been defined in (|12|) and 
where we have substituted (jblTf) . In equihbrium it is straightforward to 
check that (|2()j) coincides with the bath temperature by integrating by parts 
twice the integral in the denominator of ()26() . 

/oo pco 
dxV {x)exp{-pV{x)) = (3 dx{V (x))^ exp{-(3V{x)) 
-oo J —oo 

(27) 

In the aging regime in a partial equilibration scenario the effective temper- 
ature p6|) is related to the energy at time through the relation (|23|) . 



5 The homogeneous potential model 

An example of the GOM where many quantities can be easily worked out is 
the case of an homogeneous potential of the type, 

V,{x) = ^x'P , (28) 

with p a positive integer. The probability distribution Q is given by, 

P{AE) = {27Tk^A^h2p)-^ exp(-^ ^) , (29) 



2p 



with /ifc = a;^^^-!). By definition E = NVp = N{k/2p)h,pj^i. Similar re- 
sults to can then be obtained for the distribution of changes P{Ahk) 
for generic observables hk which lead to a hierarchy of coupled dynamical 
equations similar to (fT(1]) . These equations can then be studied using gen- 
erating functional techniques similar to those developed in other solvable 
spin-glass models |22j. Only for the harmonic case p = 1 the equation for 
the energy ()1U|) is closed and Markovian as its time evolution only depends 
on the energy. 

The expression for the effective temperature ()26() for the homogeneous 
potential model is given by, 

^^=^^'^^-(2^-1)^^ ■ ^ ^ 

In the partial equilibration scenario all observables are functions of the en- 
ergy E of the hypersurface over which the system partially equilibrates. The 
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relation between the value of T^^{tuj) and the energy E{tw) can be easily 
derived from ()23() . Introducing ()28() in @ we get, 

Zi = 1" exp(-/3Fp(x)) = (-^)^|_^dyexp(-y2p) , (31) 

yielding the free energy F = const — A^riog(^i) = const — {NT /2p) log(T) 
and the following expressions for the energy and entropy, 

_ NT 

^-^-^ ' ^^^^ 

5 = -^ = ^ + ^^^ = ^^M^ + const , (33) 

dT 2p 2p 2p ^ \ J 

giving, 

1 = (^) = ^ (34) 

T^^itw) ^ dE )E=E{t^) 2pE{t^) ■ ^ ' 

The case p = 1 corresponds to the harmonic oscillator and gives the well 
known equipartition relation T^g^^t^)) = 2E{tw)- As remarked in the para- 
graph following (|14|1 the dynamical evolution of the energy can be solved in 
general by knowing the relation 1)34(1 between the energy e = E/N and the 
effective temperature T^^ {tw) = A(i^). In this case it is possible to derive 
the following asymptotic behavior for the energy, 

e(t) ^ — TTT- + subleading logarithmic corrections . (35) 
^ ' (log(t))f ^ ^ ' 

The effective temperatures (jnU)) can be also derived using the fluctuation- 
dissipation relation (|24|) or a set of microcanonical relations describing ob- 
servable changes. The possibility to obtain the same value of the effective 
temperature by using three different approaches (the FDR H24|) , the FT (|25() 
and the microcanonical rates for observables) has been explicitly shown for 
the harmonic case p = 1 |S] . In Ref . |S] the main assumption was the validity 
of the partial equilibration scenario. Because the partial equilibration sce- 
nario also holds for the GOM, the main conclusions of Ref. |5j are expected 
to hold also for the present more general case. 

The procedure to derive the FDR (|24|) entails the computation of the cor- 
relation and response functions for the magnetization M defined as M{t) = 

I N 

C{t,ty,) = —'^Xi{t)Xi{ty,) , (36) 

i = l 



where ^(t^,) is an impulse field coupled to the magnetization M at time t^. 
At T = the response (|37() is finite due to the shift of the energy levels 
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induced by the field. For the harmonic case p = 1, dynamics is closed and a 
simple expression can be derived for the effective temperature ()24|) . 

T™«(M.) = 2e(U + -^^ , (38) 

where f{tyj) is a function that asymptotically decays as Two re- 

markable facts emerge from 1) Ea.(|H5|) only depends on the age 

at all times t > t^, therefore T™^(i,t^) = T™^(ii„) characterizes the 
aging state of the system at time t^', 2) The second term in the r.h.s of 
(jSHJ is subdominant respect to the first term. Using ()34() this leads to 
Tg™^(f^) — > T^^{tu)) — > 2e(s) so both the effective temperature derived 
from the FDR and the FT coincide. 

The obtain the effective temperature it is sometimes useful to construct 
the so-caUed fluctuation-dissipation (FD) plots US HI Hi |H|. The FD 
plots for the homogeneous model can be worked out as follows (the same 
construction holds for the GOM). As C{tw,tw) = x^{tw) = ^2(i«)) is time 
dependent then it is convenient to normalize the correlation C{t,tw) by the 
autocorrelation value taken at the lowest time C(t,tw) = C{t,tw)/C{tw,tw) 
and plotting the integrated-response xii^tw) = It^ dsR{t, s) as a function 
of C for fixed and varying t. The resulting asymptotic curve is then 
expected to have a the form of a straight line, x(C') = Xo(l ~ C) where 

Xo is the equilibrium susceptibility at zero temperature xo oc /3 p . These 
straight FD plots are characteristic of the one-step behavior observed in 
structural glasses. Fig. |31 shows the resulting FD plot for the p = I case as 
obtained from numerical simulations of the model. 



6 The wedge potential model 

An interesting example of the GOM is the case where the first derivative 
of the potential V'{x) is not continuous at the ground state configuration 
X = 0. In this case the expressions H6l7p need to be reconsidered. The 
classical example for a potential of this type is the wedge potential model 
defined by, 

V{x) = k\x\ . (39) 

This potential is depicted in Fig. 01 The statics for this model is straight- 
forward and the partition function is given by Zi = 2 /{(3k). The internal 
energy and entropy are given by the expressions (|32I33)) with p = 1/2. At 
equilibrium the one-oscillator probability density q'^°'{x) is given by, 

g<=q(x) = ^exp(-/3fc|x|) . (40) 
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A 

C 



Figure 3: FD plot for the harmonic model {p = 1) obtained from Monte 
Carlo simulations at T = with N = 1000 oscillators, k = A = 1, magnetic 
field intensity equal to 0.01 and = 1, 10, 100 (circles, triangles up and 
stars respectively). It represents xiC) and the straight line is the theoretical 
asymptotic prediction for large values of ty^. Data were averaged over 200 
dynamical histories. 




-1 -0.5 0.5 1 



X 

Figure 4: The wedge potential model (|39|) with k = I. It can be seen as a 
special case of the homogeneous model of Sec. Elwith p = 1/2. 
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To solve the off-equilibrium dynamics of this model we proceed similarly 
as was done in Sec. for a general function V{x). The main difference now 
is that the potential (|39|) is not differentiable at x = 0. While the expression 
for is well defined {{V'{x)y = k"^ is continuous for all x) the expression 
for Mae is not as the second derivative V"{x) is discontinuous at x = 0. 
The expression however can be guessed by noticing that V"{x) = a6{x) and 
using the relation, 

dxV"{x) = a = V'{oo) - V'{-oo) = 2k . (41) 

From (|6|7j) this gives, 

Mae = kA'^Ji^ = kA'^qtiO) , (42) 
where qt{x) is the one-oscillator probability density, 

1 ^ 

9*(^) = ^E'5(^-^*W) • (43) 
1=1 

The result (|42|) together with the relation, 

alj, = k'A^ , (44) 
gives the final probability distribution, 

{AE - kA'^qt{0)f' 



P(AE) = (27r/c^A2)"^ exp 



2A;2A2 



(45) 



This result can be alternatively derived doing a more elaborated but better 
controlled calculation that we do not consider interesting enough to repro- 
duce here in detail. The effective temperature (^1)) can be obtained using 
the FT p5l26|l. 

Using the thermodynamic relation (|32|) with p = 1/2 we get e = E/N = 
k\x\ = T. If we now assume that ()4Up holds in the aging regime by replacing 
P with Pesitw) = ^/Tcs{tw) then, in a partial equilibration scenario we get, 

qtJx) = ^^^^eM-Pesit^)k\x\) , (47) 

which gives (0) = (A;/5eff(i«)))/2 = k/{2e{tw))- Substituting this last result 
in H46|) we obtain, 

T!^itw) = = e{t^) , (48) 
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which coincides with the result derived using the thermodynamic relations 
H23|) and 1)33(1 with p = 1/2. Finally, we mention that the asymptotic de- 
cay for the energy derived for the homogeneous potential model (|35|) also 
holds in the wedge model substituting p = 1/2, i.e. e{t) ~ l/\/log(t). In 
Figs. I5|6|7I we show some numerical results for the wedge model obtained 
by doing Monte Carlo calculations. These have been done at zero temper- 
ature for N = 1000 oscillators with parameters A; = A = 1 and starting 
from a random initial configuration with initial coordinates Xi chosen from 
a Gaussian distribution of zero mean and unit variance. These simulations 
are useful to verify the main predictions. Fig. El shows the time decay of 
the energy for the wedge model compared to the harmonic model. Fig. [S] is 
a test of the main result ()45|) . The energy change distributions have been 
evaluated at three different waiting times = 10, 100, 1000 and fitted to 
a Gaussian distribution of variance fc^A^ = 1 where the mean qt^{0) is the 
fitting parameter. This value is then compared with the value of the energy 
e{tw) to check the identity (|l8|) . The results are shown in table I. Fig. Q 
shows the one-oscillator probability density in the aging regime qt^{x) at 
three different waiting times compared to the expected result (|47() . 

7 Conclusions 

In this paper we have introduced a new family of exactly solvable models 
characterized by zero-temperature relaxation determined by entropy barri- 
ers and partial equilibration. Generalized oscillator models (GOMs) offer a 
conceptual framework to develop a statistical description of the aging state. 
The interesting property of this class of models stems from the validity of the 
partial equilibration scenario [5. , : dynamics is ergodic and microscopically re- 
versible when the system is constrained to move along the constant energy 
hypersurface. We have then computed the probability distribution of en- 
ergy changes P(AE) that characterizes the spontaneous relaxation process 
at zero temperature. The spontaneous process is not thermally activated 









10 


0.335 


0.28 


100 


0.180 


0.178 


1000 


0.132 


0.136 



Table 1: Effective temperature in the wedge model (|48|) obtained in two 
independent ways. Second column: From the value of the energy e{tyj) 
shown in Fig. El Third column: From the value of qt^ (0) obtained by fitting 
the expression H45|) to the numerical distribution P{IS.E) shown in Fig. El 
Both values asymptotically coincide for large values of tw 
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log(t) 

Figure 5: Monte Carlo studies of the harmonic model (continuous line) and 
the wedge model (dashed line). Simulations have been done at T = with 
N = 1000 oscillators and A: = A = 1. Main: Energy decay e(t) in both 
models. Inset: e(t) multiplied by (log(t))^/^ (harmonic model with p = 1, 
wedge model with p = 1/2) as a function of log(t). In the large t limit both 
curves converge to a constant value. 

but determined by the fact that the system has been prepared in a non- 
equilibrium aging state. Using a fluctuation theorem for the aging state it 
is then possible to derive analytic expressions for the effective temperature 
(|26|) without the need to solve the dynamical equations for correlations and 
responses (jUJ). The quantitative description of the spontaneous process in 
terms of a fluctuation theorem valid in the aging state has been recently 
proposed to be at the root of the intermittency phenomenon observed in 
glasses and colloids |2l IH] . 

Two classes of models have been studied in detail. The homogeneous po- 
tential model in Sec. [Sand the wedge potential model in Sec. El Particularly 
interesting is the latter where the first derivative of the potential is singular 
at the ground state configuration. In this case, the effective temperature 
depends on the value of the one-oscillator probability distribution at 
the value of the singularity of the potential. The present studies can be 
extended to other interesting potential functions and some results will be 
presented in the future. 

A salient feature of the GOM is the Gaussian shape of the distribution 
P(AE) in large N limit. This is a direct consequence of the dynamics of the 
model © which is of the mean-field type as oscillator correlations do not 
enter the analytical expression of P{AE). In this regard, the validity of the 
partial equilibration scenario in the GOM is consequence of the mean-field 
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Figure 6: Probability distribution P{AE) (|45|) for the wedge model numer- 
ically evaluated at T = 0, iV = 1000, A; = A = 1 for = 10, 100, 1000 
(circles, triangles up and stars respectively). The continuous lines are the 
fitted Gaussians with Ma_e = Qt^i^) &s fitting parameter. The effective 
temperature ()46() is then compared with that obtained from the energy in 
Table I. 




X 



Figure 7: One-oscillator probability distribution qt^{x). The same parame- 
ters and symbols as in Fig. El Note the presence of the exponential tails in 
the distribution as expected from (|47)) . 
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character of the dynamics. The extension of these ideas to non mean-field 
dynamical models is an open problem. Nevertheless, despite of the fact that 
the present ideas have been derived from the study of mean-field systems, 
we do not foresee conceptual limitations in their adaptation to spatially cor- 
related dynamics. A conceptual description of the aging state in terms of 
heat exchange processes jU E] could be achieved in terms of a spatially fluc- 
tuating effective temperature that would describe local fluctuations in the 
rate of energy relaxation in the system. At difference with bath tempera- 
tures, effective temperatures have to be considered as fluctuating intensive 
variables as they describe energy exchange processes (either heat releasing 
or work releasing -through mechanical stresses-) occurring over nanoscale 
spatial regions where energies are not macroscopic but of the order of few 
ksT. Would be very interesting to explore the possible connection between 
these ideas and the existence of spatial heterogeneities that have received 
considerable attention during the past years (see for instance [211 HH ) ■ Es- 
tablishing a thermodynamic description of these heterogeneous excitations 
is probably an important step towards their understanding. 
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